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Abstract 



It is argued that it is valid to use QCD sum rules to determine the scalar and 
pseudoscalar two-point functions at zero momentum, which in turn determine 
the ratio of the strange to non-strange quark condensates Rsu = with 
(g = u,d). This is done in the framework of a new set of QCD Finite 
Energy Sum Rules (FESR) that involve as integration kernel a second degree 
polynomial, tuned to reduce considerably the systematic uncertainties in the 
hadronic spectral functions. As a result, the parameters limiting the precision 
of this determination are Aqcd, and to a major extent the strange quark 
mass. From the positivity of Rsu there follows an upper bound on the latter: 
m7(2 GeV) < 121 (105) MeV, for Aqcd = 330 (420) MeV . 
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1 Introduction 



The ratio of strange to light-quark vacuum condensates, Rsu = is a 

key QCD parameter measuring flavour SU(3) symmetry breaking in the vac- 
uum [1]. It is also an important quantity that enters in many QCD sum 
rule applications, e.g. baryon mass determinations, the Goldberger-Treiman 
discrepancy in SU{3) x SU{3), etc. [2j-[3j. In addition, this ratio is related 
to two low energy constants of chiral perturbation theory [4j, which in turn 
determine the next-to-leading order corrections to the Gell-Mann, Oakes, 
Renner (GMOR) relation. As a result of this importance, many attempts 
have been made in the past to determine the numerical value of this ratio, as 
well as to improve its accuracy [2], |5]-|6]. Improvements in the QCD sector 
have been possible due to state of the art results for the relevant two-point 
functions at higher order in perturbation theory, as well as to a better un- 
derstanding of how to deal with logarithmic quark-mass singularities. Better 
accuracy in the strange quark mass and in Aqcd is still required. A seri- 
ous limiting factor, though, has always been the lack of direct experimental 
information on the hadronic spectral functions entering the QCD sum rules 
used to extract Rgu- While data on hadronic r-lepton decays has allowed for 
a simultaneous determination of the (light) vector and axial-vector spectral 
functions, this is not yet possible for the scalar and pseudoscalar counter- 
parts which determine Rgu- Even if all scalar and pseudoscalar resonances 
were to be firmly established, a reconstruction of the hadronic spectral func- 
tion would remain model-dependent to a large extent. In fact, inelasticity 
and non-resonant background are hard to model correctly. 
In this paper we argue that it is valid to use QCD sum rules to determine 
the scalar and pseudoscalar two-point functions at zero momentum, which 
in turn determine the ratio Rsu- These sum rules actually fix the differ- 
ence between the true ip{5){0) and its perturbative piece. In an attempt 
to reduce systematic uncertainties from the hadronic sector we introduce a 
new set of Finite Energy QCD sum rules (FESR) to estimate the scalar and 
pseudoscalar two-point functions at zero momentum . These FESR involve 
as integration kernel a second degree polynomial with two free parameters. 
These are determined by requiring the vanishing of the spectral function at 
the position of the first two resonances in each channel. As a result of this, 
the numerical importance of the hadronic contribution to the FESR is con- 
siderably reduced. In fact, it becomes roughly an order of magnitude smaller 
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than the QCD counterpart. In addition, the latter turns out to be dominated 
by the purely perturbative QCD (PQCD) piece; the higher order in terms 
as well as the condensates add up to a negligible contribution as a result of 
partial cancellations. The results show a very good stability against changes 
in the upper limit of integration over a wide range of energies. Sensitivity to 
Aqcd, and most particularly to remains somewhat high, and becomes the 
limiting factor in the accuracy that can be achieved for Rsu- This is roughly 
at the 20 % level. Nevertheless, with the hadronic uncertainties well under 
control in this approach, future reduction of the errors in and Aqcd will 
allow for a more accurate determination of Rsu, almost free from hadronic 
systematic uncertainties. 

The high sensitivity of the results for 4'(5){0) and Rgu to the strange quark 
mass are used to derive an upper bound for this quantity from the require- 
ment Rsu > 0. Finally, the unphysical low energy constant of Chiral Pertur- 
bation Theory related to Rgu is estimated. 



2 Low energy theorem 

We introduce the correlator of vector and axial-vector divergences 

^(5)(g')ll =^ I d'x e^"^ < \Tid'^Ux)\l , d"^ 4(0)11)1 > , (1) 

where d^Jf^{x)\l = {nij T ^i) '■ Qji^) i {I5) Qi{x) '■ is the divergence of the 
vector (axial- vector) current. At zero momentum, a Ward identity relates 
the subtraction constants 'ip(5){0)\l to the quark condensates [1], [Zj-[H], viz. 

A^mi = -K- T m,) (ft-V', T . (2) 

In the determination of < > we shall be using i = u,d and j = s, as well as 
the approximations rris » rriu^d, < uu >^< dd >. From the time-ordered 
product in Eq.(l) and using Wick's theorem one would get normal-ordered 
operators in the low-energy theorem Eq.(2). However there are mass-singular 
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quartic terms in perturbative theory as well as tadpole contributions. In fact, 
to lowest order in the MS scheme 



^"m = -i^-t(l + ln^), (3) 

(0 |m,.-(0).(0)| 0)^ = -^m^(l + In 4) • (4) 

If the quark condensates in Eq. (2) would be considered as minimally sub- 
tracted instead of normal-ordered, then the perturbative quartic mass cor- 
rections would cancel [8] -[11], the low-energy theorem would make sense and 
the simple functional form of Eq.(2) would follow. All of this still holds after 
introducing gluonic corrections [5]. 

If the condensate is to be calculated from QCD sum rules, there is a subtle 
point concerning the renormalization of the operators in Eq.(l), first pointed 
out clearly in [B], which we discuss in the following. In QCD the correlator 
?/'(5)(g^) is of the general form 

where A and B are constants related to external renormalization, and ip{^){q'^) 
is the two-point function without the first order polynomial, which has been 
factored out. In the MS scheme of QCD perturbation theory, for instance, 
and at lowest non-trivial order in the strong coupling constant, the correlator 
at large —q^ is given by 



— - 34L + QL^ - 24C(3) 
L) + 0{asml,ml) , (6) 



2 2^L /lO RT , "^(g ) 
-g2^oo lOTT^ I Vr 



+ 0{al{q')) 



4 12 

lOTT^ 



where L — In and ({z) is the Riemann zeta-function. The correlator at 

zero momentum can be formally written in terms of a QCD Finite Energy 
Sum Rule (FESR), which follows from Cauchy's theorem in the complex 
energy (squared) plane (see Fig.l), i.e. 

where Sth is the hadronic threshold (e.g. M]^), and the contour integral is 
performed over a large circle where the exact ip(5){s) can be safely replaced by 
its QCD counterpart t/'^^^(s). To leading order in chiral-symmetry breaking, 
i.e. to order 0{ml), the constant A in Eq.(5) vanishes, and the hnear term 
in does not contribute to the integral in Eq.(7). 





lm(s) 
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Re(s) 









Figure 1: Integration contour in the complex s-plane. 

This means that in this case ■?/'(5)(0), as well as the non-normal ordered con- 
densate in the low energy theorem, Eq.(2), can be determined unambiguously 
from a FESR. At the next order, i.e. keeping terms of order O (m^), and 
using Eq.(5) in Eq.(7) leads to 
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V^5(0) = -- Irmp^is) + i - i^r'is) (8) 
Jsth s 2m J s 

C{\so\) 

where ijjf'^^s) is defined in Eq.(5) . The sum rule then relates the non- 
normal ordered condensate to an integral over the hadronic spectral function 
and a contour integral over the non-polynomial part of the QCD correlator 
ijj^'^^{s). Whereas the product of quark mass times normal-ordered conden- 
sate is a renormalization invariant quantity, this is, however, no longer true 
for the non-normal-ordered condensate. Only at leading order in is the 
quark condensate directly related to a physical quantity. In full QCD, how- 
ever, its value, just like the QCD coupling and the quark masses, depends 
on the renormalization scale and on the renormalization scheme employed. 

Numerically, the quartic mass terms are potentially relevant only for the 
determination of the strange quark mass. In fact, we find a-posteriori that 
the corrections of order to ip{5){0) are at the level of only (1 — 2)% of the 
leading terms. Hence, the subtleties of renormalization discussed above are 
largely academic. 



3 The ratio Rsu 

One possible way of determining the ratio Rgu is to use the auxiliary ratio 



P ^ ^5(0)1^ _ 1 rris + mu , < ss > . . 



where < uu > ~< dd > will be assumed in the sequel. In fact, if the 
subtraction constants and the quark masses are determined independently, 
e.g. from QCD sum rules, then Rsu follows. Using current values of the 
quark masses [12] gives 
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< ee > 

Rsu = — (10) 

<uu> 

Since Raa is expected from current algebra to be of order 0(10), this method 
would result in a very large uncertainty in Rgu unless the subtraction con- 
stants were to be determined with extreme accuracy. Due to this, an alter- 
native procedure, first proposed in pJij, consists in using instead the ratio 



R - nil 



which leads to 



(12) 

{uu) 1 - RvA 

This method was used in [5] to obtain both subtraction constants from 
Laplace transform QCD sum rules to four loops with the results 



^5(0)1^ = (3.35 ± 0.25) X 10"^ GeV^ , (13) 

7/,(0)|^ = -(1.06 ± 0.21) X 10"^ GeV^ , (14) 

Rsu= ^''^ =0.5 ±0.1 , (15) 

<uu> 

and the following value of the invariant strange-quark mass 

fhs = 140 ± 10 MeV , (16) 
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for Aqcd in the range Aqcd — 300 — 350 MeV. Not included in the above er- 
rors are the uncertainties due to hadronic spectral function modelling, which 
could be large jT] 



The connection between the quark condensate ratio Rg^ determined from 
QCD sum rules, e.g. through Eq.(12), and the one entering chiral pertur- 
bation theory has been discussed in [6]. In the framework of the latter, Rgu 
depends on an unphysical low-energy constant H2 through the relation 



Rsu = — = 1 + 3/i^ - 2fiK - /i„ + T^(M| - M-^){2Ll + H;) , (17) 

<qq> K 



where < gg > is the average of the up- and down quark condensates, Lg is 
a (physical) low-energy constant in the chiral Lagrangian to next-to-leading 
order [T3, and 



Ml Ml 

with the chiral renormalization scale. The constant Lg at a scale equal 
to the rho-meson mass has been estimated in chiral perturbation theory to 
next-to-leading order with the result [B] 



Ll{v^ = Mp) = (0.88 ± 0.24) x 10"^ , (19) 

while a determination at order 0{p^) gives [16J Ll(u^ = Mp) = (0.62±0.20) x 
10"'^. The unphysical low energy constant H2 has been estimated in ||6j as 



H;{iy^ = Mp) = -(3.4 ± 1.5) x 10"^ . (20) 
Both low energy constants determine the size of the next-to-leading order 
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chiral corrections to the GMOR relations in SU{2) x SU{2) and SU{3) x 
SU{3), 6t, and 6k respectively, defined as 



{rriu + rud) < uu + dd >= -2/^M^(l - 6^) , (21) 



nis < ss > {1 + = -2flM],{l - 5k) , (22) 



where the physical values of the pseudoscalar decay constants are = 92. 4± 
0.26 MeV, and fK/fn = 1-22 ± 0.01 [12]. To next-to-leading order one has 

nsi 



Ml 

6^ = 4-f{2Ll~H;) and 6k = j^S^. (23) 



4 Finite Energy QCD Sum Rules 



We consider first the pseudoscalar correlator, ip^lq'^), which exhibits a pole 
and a cut in the complex energy (squared) plane. Cuachy's theorem reads 

^5(0) = 2/iMi + ^ £ ^ V^5(^) , (24) 

where the closed contour C comprises the cut across the real axis and the 
circle of radius |so| (see Fig. 1). Introducing an integration kernel of the 
form 



A5(s) = 1 - ao5 s - ai5 , (25) 
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where aos, and 015 are free parameters, the two-point function at zero mo- 
mentum becomes 



1 Z""*" r/s 

ij.iO) = 2/|.Mi A5(M|) + - / -A,{s)Im^,{s)\RES 
X f ds 

+ — f — A,{s) ^^5){s)\qcd . (26) 

^71"^ Jc(\so\) "S 



The free parameters ao5, and ais will be chosen in such a way that A5(Mf) = 
A5 (M|) = 0, where Mi,2 are the masses of the two radial excitations of the 
kaon. This procedure will reduce considerably the numerical importance of 
the resonance contribution to 1^5(0), thus reducing the systematic uncertain- 
ties that plague the hadronic sector. For the scalar two-point function at 
zero momentum one finds 



^(0) 



1 /"'"^ ds 1 r ds 

= - — A{s) Im ^{s)\res + 7r- f - ^(^) ^i^)\QCD , (27) 



where A(s) is a second degree polynomial as in Eq. (25), and it will also 
be constrained to vanish at the position of the two resonances in the scalar 
channel. 



The two-point function V'(5)(Q'^) has been known in PQCD up to four-loops 
for quite some time ^7\. Recently, the PQCD second derivative of ip(5){q^) to 
five loops has been computed in [TS]. Integrating this result twice gives the 
five-loop expression for the two-point function up to polynomial terms . The 
latter do not contribute to the integrals around the circle in the s-plane. The 
remaining terms in the QCD expression for ip{5){q'^), i.e. the higher orders 
in fris and the quark and vacuum condensate contributions, may be found in 
To compute the QCD contribution we define 



^{5){so)\qCD = 



1 



ds 

C{\so\) S 
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A(5){s) ij(5){s)\QCD 



(28) 



where V'5(s)|qcd = iK''^)\QCD for the purely gluonic piece, but 4(so) 7^ H^o) 
on account of A5(s) ^ (095 7^ ^-nd 045 7^ 04). Using the expression 

for ip{5){q'^) to five-loop order in Eq. (28) and performing the integration 
around the circle in the complex s-plane we find the following purely gluonic 
results 



\5)\lLOOP 



167r2 



m^(so) [6 So - 3 ao(5) So - 2 ai(5) So] , (29) 



\5)\2LOOP — 



167r2 



m^(so) 



TT 



38 

46 So - 20 ao(5) «o " y «i{5) 4 



(30) 



f I l — 2/ \ r'^s('5o)-|2 I 

\5)hLOOP = -Y^Wl,(So)[ ^ J <;So 



, TT 



+ 190-51(— -1) 



«0(5) 



ai(5)- 



2 

9631 



9631 
9631 



- 105 C(3) 

- 105 C(3) + 95 



190 

24 1<'5«3) + ^ 



- 51 (^-i) 

^ 6 9' 



(31) 



>(5)|4LOOP — 



167r2 



— 2/ \ r'-^'*('^o) 13 I 



TT 



.4. + ^2i'^ - f C(3)) 



(1374 + f)(^!-l) 



00(5) - 



Ai + 6 



,4781 475 



18 



X C(3))-(1374 + f)(^4, 



ai(5)- 



A + 4( 



4781 
1^ 



g sw^ V 2 ^ ^ 6 9^ 



(32) 
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6. 



{5)\5LOOP 



1 — 2/ \ rQis(So)i4 

— m,(so) f 1 



TT 



■So 



167r 

X (_-l) + 120i/,(--- + l, 



^0(5) 



vr 



2 

2 1 



ifi — H2 



(6H3 - 12HA (—--) + 120 {- + — 

V 3 4; V g ^) 5v^20 24 16' 



0-1(5)- 



/^i - - (QHs - 8H,) + 120 H, 

6 9 



X 



120 54 ^ 8r 



(33) 



where for three quark flavours Ai = 2795.0778, Hi = 33532.30, H2 = 
-15230.645, H^ = 3962.4549, H^ = -534.05208, and H5 = 24.171875. The 
constants Hi enter in the expression of the two-point functions to five loops 
as 



^(5)(g^)|5LOOP = Y.H^L\ (34) 

167r- TT ^ — ' 

i=l 

where L = ln(— g^//i^), and the above expansion is up to (unknown) terms 
not muhiplying logarithms, which do not contribute to 5(5). The remaining 
QCD contributions to 5(5) (higher order in and vacuum condensates) have 
also been calculated, but their total contribution is at the level of 1 — 2% of 
the sum of the gluonic terms in the wide range Sq — 2 — 6 GeV^. 



Turning to the hadronic sector, the spectral function in the pseudoscalar 
channel, Imip^{s)\HAD involves the kaon pole, plus two radial excitations, 
the K(1460) and K(1830) both with widths of about 250 MeV. We follow the 
procedure outlined in [19], where the resonance part of the spectral function 
is written as a linear combination of two Breit-Wigner forms normalized at 
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threshold according to chiral perturbation theory. The latter incorporates 
the resonant sub-channel K*{892) — tt which is important due to the narrow 
width of the K*{892). This gives 



1 ds 

S5iso)\HAD = 2flMl /\^{Ml) + - — A5(s) /mV5(s)|iJ£5 

Jo S 

= S5{so)\poLE + S5{so)\res ■ (35) 

For the scalar channel there is experimental data on Kn phase shifts [20] 
that can be used to reconstruct the spectral function 



-Im^{s) = -^ /^'~'^^^'-'-^ \dis)f , (36) 

TT JZTT^ S 

where s± = [M^ ± M^)^, and d{s) is the scalar form factor. One can use 
the method of [141, based on the Omnes representation, to relate d{s) to 
the experimental phase shifts. A posteriori, the numerical importance of 
the resonance contribution to 5(5)(so)|//yiz) is one order of magnitude smaller 
than the gluonic contributions on account of the integration kernel A(5)(s). 
Hence, a simpler parametrization in terms of two Breit-Wigner forms, prop- 
erly normalized at threshold with |(i(s+)| ~ 0.3 GeV^, is equally acceptable. 
We thus include the i^^^(1430) and the Kq{1950) with masses and widths 
Ml = 1.4 GeV, Ti = 290±21 MeV, and M2 = 1.94 GeV, T2 = 201±86 MeV. 
The function 6{sq) in this channel can then be written as 

1 r° ds 

S{so)\had = - — A(s) Im iP{s)\res = S{so)\bes ■ (37) 
Jo s 

Requiring A(5)(s) to vanish at resonance determines 00(5) and 01(5) with the 
result 
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ao = 0.777 GeV-^ ai - -0.136 GeV"^ 
ao5 = 0.768 GeV"^ ais = -0.140 GeV^^ 



(38) 



The values of these coefficients are very similar on account of the similarity 
between the scalar and pseudoscalar resonance masses. 



5 Results 



In order to compute the QCD contribution to the scalar and pseudoscalar 
two-point functions at zero momentum we need as input the invariant strange- 
quark mass rhg and the QCD scale Aqcd, which are strongly correlated. To 
obtain the running quark mass and strong coupling constant it is only nec- 
essary to use the four-loop expressions, which for three quark flavours are 



1 + (290 - 256LL) 



1 1 
729L 



+ 



550435 



80 



C(3) 



1062882 729 
- (388736LL - 106496LL^) 

126940037 256 



531441 



1 

Z2 



+ 



+ 



1162261467 177147 
611418176 112640 



128 7520 

Pi + 74 + C (3) 

19683 ' 531441 ^ ^ ^ 

335011840 



387420489 531441 



C(3) LL 



387420489 



-LL' 



149946368 
1162261467' 



LL' 



1 / 1 



(39) 
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TT TT \ TT 





^(In'L-lnL- 1) + 



+ 0{l/L'), (40) 



where 



(41) 



TT (3iL ' 

L = ln(so/A|cD), LL = InL, /3i = -9/2, /^s = -8, p3 = -3863/192, and 



281198 890 

/?4 = C 3 , 42) 

4608 32 ^ ^ ^ ' ^ ^ 

with 74 = 88.5258 [21], and mj is the invariant quark-mass. The terms 
of order 0(^j^) above are known up to a constant not determined by the 
renormahzation group. However, we have checked that our final results are 
essentially insensitive to the inclusion of terms of this order in ag and mj. 



Since we are dealing with three quark fiavours, it is simpler to determine 
^QCD from the strong coupling obtained from r-decay [12], [22]: q;s(M^) = 
0.31 — 0.36, which gives Aqcd = 330 — 420 MeV. Recent determinations 
of the strange quark mass from various QCD sum rules [5], [18], [23] give 
values in the range m7(2 GeV) ~ 80 — 130 MeV, which translates into 
rfis ^ 100 — 170 MeV after using the above values of Aqcd- The two-point 
functions at zero momentum are given by 



m\: = S\REs{so)+6\QCD{so) , (43) 
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^5{0)\l = 4|polb(so) + S5\res{so) + S5\qcd{so) , (44) 

where the various S's above are computed from Eqs.(29)-(33),(35) and (37); 
as mentioned earher, the QCD contributions to 5(5) (sq) from higher orders in 
TTis and from the vacuum condensates add up (due to partial canceUations) 
to less than 1% of the sum of the gluonic terms. In Fig.l we show the results 
for i/j^lso) (curve (a)), and "0(0) (curve b) as a function of Sq for the reference 
value of the invariant quark mass nis = 100 MeV (m7(2 GeV) ~ 80 MeV), 
and Aqcd = 330 MeV. As seen from this figure the results are fairly stable 
in the wide region sq — 2 — 6 GeV^; a similar stability is obtained for Aqcd = 
420 MeV. For values in the range Aqcd = 330 - 420 MeV we find 

^5(0) = (0.39 ± 0.03) X 10"^ GeV^ (m, = 100 MeV) , (45) 



^(0) = -(0.95 ± 0.25) X 10"^ GeV^ (m, = 100 MeV) , (46) 
which using Eqs.(ll)-(12) leads to 



Rsu = i=rv = 0.6 ± 0.1 (rfis = 100 MeV) . (47) 
[uu) 

The uncertainties above are entirely due to the uncertainty in Aqcd, as 
has been kept fixed at the indicated reference value. The ratio Rgu exhibits 
a stronger sensitivity to the value of the strange quark mass, as this enters 
in the PQCD expression of -0(5) as an overall multiplicative factor rhs^- 
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Figure 2: The two-point functions at zero momentum, ^5(0), curve (a), 
and ■0(0), curve (b), as a function of Sq, for ffis — 100 MeV, and h-qco — 
330 MeV. 



Given the similarity between the integration kernels in the scalar and pseu- 
doscalar channels, A(s) ~ A5(s), or 5\qcd{sq) ~ ^5|qcd('So), it is possible 
to obtain an approximate expression for Rsu as a function of as follows. 
From Eqs.(12) and (43)- (44) one has 
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POLE + ^5 


RES + S^IqCD + S\rES + S\qCD 


^5 


\pOLE + ^5 


\res — Sres + S5\qcd — S\qcd 





(^5 POLB 


+ 
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+ S\res 


+ 


^5 QCD + S\qCD 




^5 POLE 


+ 


^5\rES 


— S\res 


^blpOLE + RES — ^\rES 



A +B{Aqcd) 



ms{MeV) 
1 100 MeV 



(48) 



where we approximated [^sIqcd — ^Iqcd] — in the denominator of the above 
ratio, A ~ 1.15 is basically constant in the wide range sq = 2 — 6 GeV^, and 
B{Aqcd) = -0.44 (-0.68) for Aqcd = 330(420) MeV, respectively. This 
formula is accurate to within 2 - 3 %, and it allows for a quick estimate of 
Rsu for other values of the invariant strange quark mass. It also gives an 
upper bound for this mass from the fact that Rgu > 0, viz. 



162 MeV (Aqcd = 330 MeV) 
rris < { (49) 
130 MeV (Aqcd = 420 MeV) . 



These bounds translate into the following bounds for the running strange 
quark mass at a scale of 2 GeV 



r 121 MeV (Aqcd = 330 MeV) 
ml(2GeV) < <^ (50) 
[ 105 MeV (Aqcd = 420 MeV) . 

These results are in line with recent determinations from QCD sum rules 
[5J, [18j, [23], as well as Lattice QCD [24j. However, in making compar- 
isons, the strong correlation between rfis and Aqcd should be kept in mind. 
In particular, older determinations giving higher values of used mostly 
AccD ^ 100 - 250 MeV. 
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Figure 3: The ratio Rsu as a function of sq, for rfis = 100 MeV, Aqcd = 
330 MeV, curve (a), and Aqcd = 420 MeV, curve (b). 



We now turn to the imphcations of these results for chiral perturbation the- 
ory, as outhned in Section 2. Inserting our resuh for Rsu, Eq.(47), in Eq. 
(17), and using Eqs. (18) and (19) gives the following prediction for the low 
energy constant H2 

H; = -(4.3 ± 1.3) X 10^^ , (51) 

where the range = 82 — 92 MeV was used, to take into account uncer- 
tainties from higher orders in the chiral expansion [6]. 
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Using this result together with Eq. (19) in Eq. (23), the next-to-leading 
order corrections to the GMOR relation become 



6^ = 0.04 ± 0.02 

(52) 

6k = 0.5 ±0.2 , 



in good agreement with [6]. 



6 Conclusions 



We have argued that it is legitimate to use QCD sum rules to determine the 
scalar and pseudoscalar two-point functions at zero momentum. These sum 
rules actually fix the difference between the true ip(5){0) and its perturba- 
tive piece. Approaches based on traditional QCD sum rules, e.g. Laplace 
transform sum rules, are affected by uncontrollable systematic uncertainties 
in the reconstruction of hadronic resonance spectral functions. To minimize 
these uncertainties we have introduced new Finite Energy QCD sum rules 
(FESR) involving an integration kernel in the form of a second degree poly- 
nomial with two free parameters. Requiring the hadronic spectral function 
to vanish at the position of the first two resonances determines these con- 
stants, and reduces the importance of this contribution to the FESR by one 
order of magnitude. This makes '?/'(5)(0) dependent mostly on the strange 
quark mass, and to a lesser extent on Aqcd- The dependence on the ra- 
dius So of the integration contour in the complex energy plane is very mild, 
with the results for ip(5){0) showing very good stability in the wide range 
So — 2 — 6 GeV^. Our results for the scalar and pseudoscalar correlators at 
zero momentum, Eqs.(45)-(46), as well as for the ratio Rsu, Eq- (47), are 
in broad agreement with most previous determinations based on traditional 
QCD sum rules (Laplace, FESR) [2], [5]-[6j. However, it should be kept in 
mind that many of the old dterminations used much lower values of Aqcd, 
and somewhat higher values of m^. In addition, they used available PQCD 
results at the time, which were limited to two-, three- or at most four-loop 
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order, many were affected by logarithmic quark-mass singularities, and by 
uncontrollable systematic hadronic uncertainties in all cases. Our result for 
Rsu is somewhat smaller than one from a recent lattice determination j25] . 
The bound obtained for rris, Eqs.(49)-(50), is in good agreement with recent 
results from QCD sum rules |23], as well as lattice QCD [23], which point 
to values smaller than in the past. Future improvement in the accuracy of 
^QCDi and particularly in that of m^, will allow for a more precise deter- 
mination of the scalar and pseudoscalar correlators at zero momentum, and 
hence of Rsu-, almost free of systematic hadronic uncertainties. This becomes 
possible due to the introduction of the integration kernel, Eq.(25), in the 
FESR. It should be mentioned that this kernel, while vanishing at the res- 
onance peaks, it does not vanish at the point s = Sq where the integration 
circle in the complex energy plane meets the real axis. It is known that in 
some applications of FESR, e.g. in tau-decay, perturbative QCD does not 
seem to hold close to the real axis; this has led to the proposal of weighted 
FESR with weight functions vanishing at s = Sq |26]. In the application 
discussed here, though, this problem does not seem to arise. In fact, the per- 
turbative expansion appears to converge very well, and the stability region 
is unusually broad, extending well above standard values. The introduction 
of an additional integration kernel vanishing aX s = sq would not seem to 
provide any additional advantages. In any case, we have confirmed this by 
an explicit calculation. The subtraction constants '?/'(5)(0) remain essentially 
the same if a weighted kernel is added. 
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